3 0 9 6 £ 0 V QV 


TECHNICAL  REPORT  ARMID-TR-77001 


INCORPORATING  AN  INTERNAL  MASS-SPRING 
SYSTEM  INTO  THE  TRAJECTORY  EQUATIONS 
FOR  A SIX-DEGREE-OF-FREEDOM  MISSILE 

EDWARD  B.  LACHER 


APRIL  1977 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

MANAGEMENT  INFORMATION  SYSTEMS 
DIRECTORATE 

DOVER.  NEW  JERSEY 


APPROVED  FOR  PUBLIC  RELEASE.  DISTRIBUTION  UNyMITED. 


The  findings  in  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position. 


DISPOSITION 

Destroy  this  report  when  no  longer  needed.  Do  not 
return  to  the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh«n  Dele  Entered) 


READ  INSTRUCTIONS 


REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 

S RECIPIENT'S  CATALOG  NUMBER 


He  PORT  NUMBER 


GOVT  ACCESSION  NO. 


RMID-TR-77001 


[Technical  Report 


■mm  ii  i ■ Incorporating  an  internal 

MASS-SPRING  SYSTEM  INTO  THE  TRAJECTORY 
EQUATIONS  FOR  A SIX-DEGREE-OF-FREEDOM 
MISSILE  » 


« PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfa) 


7.  AUTHORfa) 


Edward  B.^Lacher  / 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  • WORK  UNIT  NUMBERS 


9 PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Management  Information  Systems  Directorate 
Mathematical  Analysis  Division 


Picatinny  Arsenal,  Dover.  NJ  07801 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

ARRADCOM 

Management  Information  Systems  Directorate 
Dover,  NJ  07801 

R MONITOI 


APRIL  M77 


IS.  SECURITY  CLASS,  (ol  thl,  report) 


G AGENCY  NAME  « AOORESS (II  dllltrm,  I fra*  Conlro/llnS  Oil  led) 


UNCLASSIFIED 


ISA.  DECLASSIFICATION/DOWNGRADING 
SCHEOULE 


IS.  DISTRIBUTION  STATEMENT  (ot  Olid  Report) 


Approved  for  public  release;  distribution  uni: 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  In  Block  20,  II  different  from  Report) 


19.  KEY  WOROS  ( Continue  on  reeeree  elde  II  neceeeery  end  Identity  by  block  ntmeber) 

Trajectory  equations 
Six-degree-of-free'dom  missile 
Mass- spring  system 
Missile  trajectory 


-This  report  describes  the  analysis  necessary  to  extend  the  existing  six-degree- 
of-freedom  missile  trajectory  equations  to  include  a missile  containing  an  internal 
mass-spring  system.  The  analysis  is  divided  into  two  main  parts,  dealfng  with  the 
internal  mas9-spring  system  alone  and  with  the  combined  projectile  and  mass  system 
This  work  may  be  considered  an  extension  of  the  basic  trajectory  analysis  documente 
in  Picatinny  Arsenal  Technical  Report  3391,  June  1966.  The  major  result  consists  of 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fR»A  Daia  BaIataW) 


HCUWITV  CLARIFICATION  OF  TMH  mOEfll«n  D—  «»»wO 

20.  (Cont'd) 

V 

'an  additional  differential  equation  describing  the  motion  of  the  internal  mass-spring 
system . to  be  added  to  the  basic  trajectory  equations . 

This  analysis  was  applied  to  extending  a six-degree-of-freedom  missile  tra- 
jectory computer  program  at  Picatinny  Arsenal , and  a brief  discussion  of  this  is 
included.  . 


‘ 


i 

( 


UNCLASSIFIED 

SECURITY  CL ASSl FIC ATlOM  Of  TMI*  D mi*  Bmlmrml) 


The  citation  in  this  report  of  the  names  of  commercial  firms  or 
commercially  available  products  or  services  does  not  constitute  official 
endorsement  or  approval  of  such  commercial  firms,  products,  or  services 
by  the  US  Government. 


HMXHOmCD 

jUSTlfWATW.. 


ftfl&edki 


8 OlSTlltt  BTI M / (tlT*  «« 

tkL  uHTut/M  wC 


□ □ 


ACKNOWLEDGMENTS 


The  author  wishes  to  thank  Bruce  Barnett  for  the  informative 
discussions  between  him  and  the  author  on  the  motion  of  the  internal  mass 
with  respect  to  the  projectile  and  on  moving  coordinate  systems  in  general . 

The  author  wishes  to  thank  Eugene  Friedman  for  his  help  in  de- 
veloping the  model  used  as  a basis  of  this  analysis. 


M 


TABLE  OF  CONTENTS 

Page  No. 


Introduction  1 

Discussion  2 

Model  2 

Internal  Mass-Spring  System  6 

Combined  Projectile  and  Mass  System  14 

Application  to  the  Trajectory  Computer  Program  18 

References  19 

Symbols  21 

Distribution  List  23 

Figures 

1 Internal  configuration  of  projectile  3 

2 Overview  of  projectile  position  8 


I 


I 

La 


INTRODUCTION 


This  report  describes  the  analysis  necessary  to  extend  the  six-degree- 
of-freedom  missile  trajectory  equations  to  include  a missile  containing  an 
internal  mass-spring  system.  It  also  describes  how  this  analysis  is  applied 
to  implementing  corresponding  modifications  to  the  most  commonly  used 
six-degree-of-freedom  missile  trajectory  computer  program  at  Picatinny 
Arsenal,  to  calculate  the  trajectories  of  such  projectiles  with  internal  mass 
systems . 

The  basic  analysis,  described  after  the  discussion  of  the  model,  falls 
naturally  into  two  main  parts . The  first  is  a treatment  of  the  internal  mass- 
spring system  alone  and  consists  of  the  derivation  of  the  additional  equation 
describing  the  motion  of  this  internal  system  only.  This  is  the  major  step 
required  to  incorporate  a treatment  of  the  internal  movable  mass  into  the 
existing  trajectory  equations.  The  second  part  deals  with  the  combined 
projectile  and  mass  system  whose  motion  is  governed  by  the  existing  basic 
trajectory  equations  as  previously  developed  in  Reference  1 . Since  the 
structure  of  those  equations  is  unchanged,  they  are  not  repeated  here. 
However,  now  that  the  total  missile  is  not  rigid,  a number  of  modifications 
to  equations  related  to  the  combined  system  must  be  made.  These  are  re- 
lated to  the  variation  of  the  center  of  gravity  and  involve  mostly  aerody- 
namic considerations. 

The  application  of  the  analysis  developed  here  to  extend  the  capabil- 
ity of  a trajectory  computer  program  to  handle  a modified  missile  with  an 
added  movable  internal  mass  has  been  completed . Since  that  computer  pro- 
gram was  based  on  the  analysis  in  Reference  1,  the  present  report  describ- 
ing analysis  for  extending  the  program  can  be  considered  an  extension  of 
that  document,  and  much  of  the  notation  is  borrowed  from  it.  Therefore, 
it  is  recommended  that  Reference  1 be  consulted  in  conjunction  with  the 
present  report. 

Although  this  work  is  geared  to  developing  analysis  to  be  applied  to 
a particular  trajectory  computer  program , most  of  the  topics  in  this  report 
consider  matters  relevant  to  any  six-degree-of-freedom  missile  trajectory 
computation . This  includes  the  final  section  in  which  a brief  summary  is 
presented  of  the  computer  program  changes  required. 
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DISCUSSION 


Model 


The  total  model  used  in  this  study  is  divided  into  a number  of  separ- 
ate parts.  The  physical  description  is  presented  first,  followed  by  the 
assumptions  made  limiting  the  scope  of  the  study.  Next,  the  coordinate 
systems  required  for  the  analysis  are  described,  and  finally  the  overall 
analytical  approach  is  summarized.  Together  these  explain  the  basic 
physical  and  analytical  preliminaries  and  provide  an  expanded  introduc- 
tion to  the  two-part  analysis  of  the  projectile  and  internal  mass  system 
which  follows. 

Physical  System 

The  system  under  consideration  in  this  analysis  is  an  axially 
symmetric  projectile  containing  a mass  which  is  constrained  to  move  along 
the  longitudinal  projectile  axis  only  and  is  connected  to  the  projectile  by 
a spring.  The  basic  physical  arrangement  and  internal  configuration  of 
the  projectile  is  shown  schematically  in  Figure  1 . The  internal  mass  is 
connected  to  the  projectile  by  a spring  (with  constant  k)  from  the  point  Q, 
fixed  in  the  projectile,  to  the  point  P,  located  at  the  center  of  gravity  of 
the  mass.  All  distances  shown  in  the  figure  are  measured  from  the  point 
Q,  positive  along  the  direction  labelled  +X„.  In  particular,  d represents 
the  displacement  of  the  internal  mass  and  d the  equilibrium  position  of 
the  spring,  both  measured  to  the  center  of  gravity  of  the  mass.  The  center 
of  gravity  labeled  "total"  includes  both  the  projectile  body  and  the  internal 
mass,  while  that  of  the  "body"  excludes  the  internal  mass. 

For  purposes  of  clarity  it  is  essential  to  differentiate  between 
the  missile  body  and  the  total  or  combined  system . When  referring  to  the 
basic  projectile  body  alone,  excluding  the  internal  movable  mass,  but  in- 
cluding any  other  components  (such  as  fuzesl  rigidly  attached  to  the  body, 
the  term  "missile  body"  will  be  used.  However,  when  discussing  the  aero- 
dynamics, the  more  customary  term  "missile  skin"  will  generally  be  used, 
but  it  should  be  taken  to  mean  the  same  as  "missile  body".  Roth  these 
terms  always  indicate  that  the  internal  mass-spring  system  is  not  included. 
When  referring  to  the  "combined  system,"  "overall  system,"  or  "total  sys- 
tem," however,  both  the  body  and  the  mass-spring  system  are  always  in- 
cluded. References  to  the  "mass-spring,"  "mass  alone,"  internal  mass," 
or  "internal"  system  indicate  the  internal  mass-spring  system  excluding 
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the  missile  body  or  skin.  The  terms  "missile"  and  "projectile"  are  used 
interchangeably  and  no  significance  should  be  attached  to  the  use  of  one  in 
place  of  the  other  . 

Assumptions 

For  purposes  of  this  study,  it  was  necessary  to  focus  only  on 
certain  effects  of  the  addition  of  the  mass-spring  system  on  the  projectile 
motion,  rather  than  to  solve  the  problem  in  complete  generality.  This  re- 
sulted in  some  simplification  of  the  task  by  making  possible  a number  of 
assumptions  Some  of  these  assumptions  were  also  necessitated  by  the  de- 
sire to  implement  the  calculations  by  modifying  an  existing  trajectory  pro- 
gram with  its  own  limitations.  These  . .mitations  either  were  considered 
acceptable  hir  the  current  study  or  the  model  was  modified.  The  assump- 
tions and  limitations  are  summarized  briefly  as  follows: 

It  is  assumed  that  the  projectile  body  is  symmetric 
about  the  longitudinal  X-axis 

It  is  further  assumed  that  the  internal  mass  is  also 
symmetric  about  this  axis  and  that  it  is  constrained  to  move  along  this 
axis  only  , with  negligible  tilting  or  rotation.  The  result  is  that  1 = I 

for  all  ol  the  systems  in  question. 

The  trajectories  are  limited  to  non-thrust  flights  only 
which  involve  no  fuel  burning  during  the  part  of  the  trajectory  to  be  ana- 
lyzed This  implies  that  no  mass  changes  or  jet  effects  need  be  considered 


The  mass  of  the  spring  is  neglected  in  the  analysis  of 
the  internal  system;  however,  one  could  account  for  it  using  methods  in- 
volving an  equivalent  internal  mass  increased  by  one-third  of  the  mass  of 
the  spring 


The  present  analysis  is  applied  to  short  trajectories  so 
that  only  the  "flat  earth"  model  ■ ' -scribed  in  Reference  1 is  needed 

Coordinate  Systems 

In  order  to  find  the  equations  of  motion  for  the  physical  sys- 
tem ir,  question,  it  is  necessary  to  establish  and  inertial  coordinate  system. 
1 . whose  origin  may  be  assumed  to  be  at  the  center  of  the  earth . A fixed 
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plane  system,  1!  (for  Horizontal  plane) . with  origin  at  the  center  of  gravity 
of  the  total  projectile-mass  combination  will  be  utilized  (as  in  Ref  1)  . The 
H-system  X-axis  lies  along  the  missile  axis,  the  Y-axis  is  constrained  to 
lie  in  a horizontal  plane,  and  the  Z-axis  is  chosen  so  as  to  form  a right 
handed  (rectangular)  coordinate  system.  Because  the  only  allowable  mo- 
tion of  the  mass  m is  along  the  projectile  axis,  the  only  non-ignorable  co- 
ordinate in  the  H-system  will  be  the  X-axis. 

Approach 

In  this  physical  system  consisting  of  two  distinct  parts  (the 
projectile  body  and  the  internal  mass)  , there  are  a number  of  possible 
approaches  to  describing  the  motions.  Since  the  object  of  this  study  is  to 
analyze  both  the  motion  of  the  internal  mass-spring  and  the  motion  of  the 
total  projectile  system,  it  is  logical  to  choose  these  as  the  two  systems  on 
which  to  concentrate.  The  first  part  of  the  analysis  treats  the  internal 
mass-spring  subsystem  alone  and  the  second  treats  the  total  system  con- 
sisting of  the  projectile  body  and  internal  mass  combined. 

In  order  to  investigate  the  motion  of  a physical  system,  it  is 
necessary  to  include  the  governing  equations  of  motion  and  the  forces  and 
moments  acting  on  the  system.  These  are  discussed  separately  for  each 
of  the  two  systems.  Because  of  the  differeing  nature  of  the  systems  and 
the  applicability  of  much  of  the  rigid-body  trajectory  analysis,  the  empha- 
sis is  different  in  each  case. 

In  analyzing  the  motion  of  the  internal  mass-spring  system, 
one  first  applies  Newton's  Second  Law  to  the  internal  mass  alone,  taking 
into  account  all  forces  acting  on  it.  The  differential  equation  derived 
through  this  process  is  the  desired  new  equation  of  motion  of  the  internal 
system.  The  forces  considered  acting  on  the  internal  mass  are  the  spring 
force,  giavity,  and  a frictional  term.  No  consideration  of  aerodynamic 
forces  (or  winds)  is  necessary  because  the  mass  is  inside  the  missile  body . 
No  moments  are  considered  because  for  the  present  analysis  the  internal 
mass  is  allowed  only  one  degree  of  freedom  inside  the  projectile,  without 
rotation.  These  considerations  are  treated  in  the  section  dealing  with  the 
internal  mass-spring  system. 

To  analyze  the  motion  of  the  total  system,  one  must  add  the 
new  differential  equation  from  the  internal  system  to  the  standard  equations 
described  in  Reference  1 . This  forms  a larger  set  of  differential  equations 
which,  when  solved  simultaneously  (numerically)  , describe  the  motions 


5 


of  both  the  internal  and  combined  systems . No  new  derivations  of  those 
standard  equations  of  motion  are  required  for  this  combined  system . 
Aerodynamic  and  gravitational  forces  and  moments  are  included,  and  no 
jet  effects  need  be  considered  because  of  the  assumption  of  no  thrust. 

The  basic  methodology  (and  computer  program)  was  used  previously  to 
calculate  the  trajectories  of  rigid  missiles  only.  The  present,  compound, 
non-rigid  system  whose  center  of  gravity  is  no  longer  fixed  in  the  system, 
requires  a number  of  modifications.  These  involve  aerodynamic  forces 
and  moments  and  moments  of  inertia,  and  are  considered  in  the  discussion 
of  the  combined  system. 

It  should  be  noted  that  it  would  entail  more  work  to  consider 
the  projectile  body  alone,  without  the  internal  mass  (instead  of  the  total 
system)  as  the  second  system  of  interest.  If  this  were  done,  the  existing 
trajectory  computer  program  would  then  need  to  be  applied  to  the  projec- 
tile body  alone,  as  if  it  were  being  flown  separately,  and  its  effects  on  the 
internal  mass  would  then  be  calculated.  This  alternate  procedure  would 
be  awkward  and  would  not  make  the  fullest  use  of  the  existing  analysis 
and  computer  program. 

Internal  Mass-Spring  System 

The  basic  equation  to^be  applied  to  the  mass-spring  system  is  sim- 
ply Newton's  Second  Law,  F = ma.  In  view  of  the  assumed  constraints 
on  the  mass  m,  the  only  direction  in  which  motion  of  this  mass  can  occur 
is  along  the  X-axis  of  the  H-system . Using  the  double  subscript  XH  to  in- 
dicate components  of  vectors  along  this  direction,  the  only  resulting  non- 
trivial equation  is 


XH 


ma 


XH 


(1) 


This  equation  must  be  applied  in  an  inertial  frame  only,  so  that  the 
term  on  the  right  must  refer  to  the  proper  component  of  the  acceleration 
of  the  mass  with  respect  to  the  inertial  coordinate  system  I.  To  make 
Equation  1 more  explicit  then,  it  can  be  re-written  as 


dJ 


XH 


= m 


I rm 


XH 


(2) 
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A few  words  of  explanation  are  in  order  to  clarify  the  notation  used 
for  the  vector  derivative  appearing  on  the  right  hand  side  of  Equation  2 . 

By  a (time)  derivative  of  a vector  relative  to  a coordinate  system  is  meant 
that  one  first  expresses  the  vector  in  components  of  that  system  and  then 
differentiates  with  respect  to  time . The  two  systems  used  here  are  the  in- 
ertial (I)  coordinate  system  necessary  to  apply  Newton's  Second  Law,  and 
the  H coordinate  system  along  whose  X-axis  the  mass  moves.  The  two  sys- 
tems are  in  motion  with  respect  to  each  other,  so  that  derivatives  of  a vec- 
tor are  in  general  different  relative  to  the  two  systems.  The  vector  deriva- 
tive in  Equation  2 is  relative  to  the  I-system,  but  only  its  XH  direction 
component  is  required . 

The  method  followed  expresses  both  sides  of  Equation  2 in  terms  of 
quantities  in  the  H-system;  that  is,  in  terms  of  the  displacement  d (along 
the  X-axis  in  this  system)  and  its  derivatives.  The  reason  for  this  approach 
is  that  one  can  most  easily  express  forces  in  the  H-system  and  that  one 
desires  an  equation  in  this  system  in  order  to  use  the  previous  analysis 
and  computer  program.  The  result  will  be  a single  second-order  ordin- 
ary differential  equation  for  d,  which  can  then  be  added  to  the  standard 
equations  of  motion  for  the  projectile  and  mass  system . 

Derivative  Term 

To  facilitate  the  conversion  of  the  right  hand  side  of  Equation  1 
to  the  H-system,  one  can  make  use  ol  Figure  2 in  which  the  origin  of  the 
H-system  is  located  at  the  overall  center  of  gravity  of  the  total  projectile. 
Thus , 


-1  -X  -w 

r = r + R (3) 

m 


expresses  the  position  vector  r of  the  center  of  gravity  (P)  of  the  internal 
mass  with  respect  to  the  I-system,  in  terms  of  the  position  vector  r of  P in 
the  H-system  and  the  vector  R indicating  the  position  of  the  H-system  ori- 
gin O All  of  these  vectors  vary  with  time  as  the  projectile  moves. 

Differentiating  Equation  3 relative  to  the  I-system, 


d.  r 
_I m 

dt 


dj  (r  ♦ R) 


V 


dj  R 


■MNtotMHI 


dt 


dt 


dt 


(4) 


The  standard  formula  for  expressing  the  time  derivative  of  a vector  A 
relative  to  a coordinate  system  I,  in  components  of  another  system  H in 
motion  with  respect  to  the  first  and  having  angular  velocity  to  with  respect 
to  it,  is  given  by 


diA  = 
dt  dt 


+ 05  x A 


This  equation  can  be  found  in  many  mechanics  texts  (among  them  Ref  2) 
and  also  in  Reference  1 . 

Applying  Equation  5 to  the  first  term  on  the  right  in 
Equation  4 yields 


dI  rm  dH  r - - dI  B 

~dT  ' ~dT  + “xr  * ~5r 


Differentiating  and  then  applying  Equation  5 to  the  term  in  square  brackets 
below , 


d?  r 
I m 


df  (dJLl-  - dlM  dI  (dHr  - A d! 
dfi  / dH **  -*  - /dHr  - A d!“ 

dF  \“d T +*xt)+  + coxr|  + — 


This  results  in 


d,  r d*  r d r dj  R 

I m H „ H . -*■  dw  -x  I 

+ 2co  x — + to  x ((o  x r)  + — — x r + — — 


Cross-Product  Terms 


To  simplify  the  cross-product  terms  in  Equation  7 , a few  ob- 
servations are  necessary . Since  the  motion  of  the  internal  mass  is  along 

the  H- system  X-axis  only,  that  is,  parallel  to  the  unit  vector  i one  is 
concerned  with  the  component  of  Equation  7 in  this  direction  only.  Using 
Figures  1 and  2 it  can  be  seen  that 

r = (d  - Dt)  I-*  (8A) 

k — * 

the  vector  r having  no  j or  components.  Similarly,  when  differentiated 
relative  to  the  H-system 


dH  r • • - 

”dT  = (d'DTMH 


dH  r _ ” r1- 

dt”  “ (d  " dt}  *h 


dH  r 

It  is  evident  from  Equations  8A  and  8B  that  both  to  x - and 

-k 

dco  -*■  — 

x r have  no  components  in  the  i direction  and  therefore  can  be 


dropped . 


Representing  the  angular  velocity  w = a>v„  i„  + a\,tI  i „ 

XH  H YH  H 


+ ooZH  k^,  it  is  clear  that 


then. 


“Xr  * W ' V “°ZH  )H  ' “YH  kH> 


to  x (03  x r)  = - (d  - Dt)  (to*  + co*H)  1 +. 


where  only  the  i component  is  of  interest. 
H 


I 


Substituting  the  results  of  these  observations  into  Equation  7 
and  multiplying  by  m yields  an  equation  without  cross  products. 


d?  r 


m 


I m 
2 


dt 


d?  R 


XH 


= m (d  - Dt)  - m (d  - DT)  (m*,H  + m^)  + m 


XH 


(10) 


D^.  Terms 


The  aim  is  to  obtain  a single  equation  for  the  single  unknown 
d in  terms  of  quantities  known  or  computed  from  the  basic  trajectory  equa- 
tions. The  variable  and  its  derivatives,  describing  the  position  of  the 
(movable)  overall  center  of  gravity  relative  to  the  projectile  fixed-point 
Q,  are  not  independent  of  d and  must  be  eliminated. 


The  origin  of  the  H-system,  0H<  lies  at  the  center  of  gravity 


of  the  overall  system . Thus  the  sum  of  the  first  moments  of  the  total  mass 
about  this  point  is  zero.  That  is, 


m (Dn  - D ) + m(d  - D ) = O 


(Hi 


Solving  for  D^, , 


D„  = 


mB  Db  + md 


m + m 

D 


(12A) 


Differentiating  twice, 


D_  = 


md 


mB  + m 


(12B) 


D„ 


md 


mB  + m 


(12C) 
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Substituting  these  into  Equation  10  yields 


di  r 

I m 


= — d + (D  - d)  

XH  mB*“  ® mE 


TV  ^yh*"‘zh) 


Derivative  of  R 

To  evaluate  the  derivative  in  the  last  term  on  the  right  in 
Equation  13,  it  is  convenient  to  call  the  velocity  vector  of  the  projectile 

- dI* 

center  of  gravity  with  respect  to  the  I system  V = — . Then  using 


Equation  5 


i (vM  = djl  M SI 

dt  \ dt  / dt  dt 


(o  xV  (14) 


Using  the  notation  V = VXH  iH  + VyH  jH  + VZH  then 


3 /-»-V 

XU  dt 


toxV  (a5YH  VZH  VYH<0ZH)iH 


Combining  these 


d''R 

dt’  XH  VXH  * (mYH  VZH  " VYH  “zH1  1161 


' 


7 


■ 1 


m. 


G„  = 


B 


B mn  + m 

13 


Substituting  Equation  16  into  Equation  13  and  using 
, the  right  hand  side  of  Equation  2 becomes  finally 


m 


d r 
I m 

dt* 


XH 


= mGBd 


+ (Dn  - d)  m G 


B 


(t°YH 


mZH] 


+ m V 


+ m (03yh  vzh 


VYH  “W 


XH 


(17) 


Forces 

Returning  to  the  left  hand  side  of  Equation  2,  it  remains 
necessary  to  find  the  forces  on  the  internal  mass  in  the  direction  of  the 
H-system  X-axis  only.  There  are  three  forces  which  will  be  considered: 
the  spring  force,  gravity,  and  a frictional  or  damping  term.  No  aerodynamic 
forces  need  be  considered  for  this  system  because  the  mass  lies  entirely 
inside  the  projectile. 

The  spring  force  is  -k(d-d  ) , and  the  gravity  term  can  be 

eq  _k 

expressed  as  the  component  of  the  gravitational  force  along  the  i direction , 

H 

which  is  denoted  by  mg  . 

XH 

The  frictional  term  is  used  for  the  purpose  of  obtaining  a 
damping  effect  on  the  motion  of  the  mass,  and  it  is  actually  a viscous 
damping  form  -bd.  Here,  b represents  the  damping  coefficient  where  the 
force  is  directed  against  the  motion  of  the  internal  mass  and  is  proportional 
to  its  velocity . Although  this  is  not  the  customary  means  of  representing 
static  or  dynamic  friction,  for  the  purposes  of  this  analysis  one  can  easily 
investigate  some  effects  of  friction  in  this  way . In  the  actual  application 
of  this  analysis,  little  use  was  made  of  this  frictional  term. 

Substituting  these  forces  into  Equation  2 using  Equation  17 
yields  the  final  result 

k (d  ~ d„J  _bd  + mi?xH  = m V + m GB  (DB  ‘ d)  lU>YH  + t°*ZH) 


eq 


+ m V , + m (o\,„  V„„  - V„„  03_., ) 


XH 


YH  ZH  YH  ZH 
(18) 
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In  the  actual  trajectory  computer  program,  the  above  equation  is  solved 
for  the  second  derivative  d,  and  this  is  reduced  to  a system  of  two  first- 
order  equations  expressing  the  derivatives  of  d and  d. 

Combined  Projectile  and  Vass  System 

The  overall  approach  to  describing  the  motion  of  the  combined  pro- 
jectile body  and  internal  mass  system  has  been  divided  into  two  parts. 

The  first  involved  the  internal  mass-spring  system  and  the  second  involves 
the  total  system.  In  each  case  the  general  equations  of  motion  are  applied 
to  obtain  the  specific  equations  desired.  Only  the  external  forces  and/or 
torques  acting  on  the  corresponding  system  are  considered.  This  proce- 
dure yields  the  required  set  of  equations  describing  the  motions  of  both 
the  internal  mass  and  the  entire  projectile. 

Since  the  basic  structure  of  the  acceleration  terms  (on  the  right  hand 

side  of  F = ma)  is  the  same  as  for  the  standard  trajectory  analysis,  no  new 
motion  equation  derivations  are  necessary.  The  resulting  equations  of 
motion  for  the  combined  system  are  simply  the  new  equation  (Eq  18)  plus 
the  original  trajectory  equations. 

In  the  preceding  sections,  the  forces  on  the  internal  mass  alone 
(along  the  direction  of  its  motion)  have  been  considered.  No  torques  were 
necessary  for  these  calculations,  because  only  longitudinal  motion  was 
permitted  with  no  rotation . It  remains  to  consider  here  the  external  forces 
and  moments  acting  on  the  total  missile  system.  Since  the  movable  mass  is 
a part  of  the  combined  system,  its  effect  on  the  combined  system  is  internal 
only.  The  only  external  forces  in  the  present  problem  are  the  same  ones 
already  treated  in  the  analysis  of  ordinary  missile  trajectories  (Ref  1)  . 
However,  some  consideration  as  to  how  the  motion  of  the  internal  mass  af- 
fects the  trajectory  of  the  overall  system  will  show  why  a number  of  modi- 
fications to  these  forces  and  moments  computed  for  the  equations  of  motion 
are  now  necessary. 

Using  the  present  approach,  one  does  not  consider  the  forces  and 
torques  applied  to  the  projectile  body  by  the  internal  mass,  because  they 
were  accounted  for  when  the  forces  exerted  on  the  internal  mass  were  con- 
sidered. In  view  of  this,  it  might  not  be  evident  how  the  present  model 
accounts  for  the  internal  mass  affecting  the  motion  of  the  combined  system . 
The  answer  is  that  as  the  mass  moves  inside  the  projectile,  it  transmits 
effects  of  its  motion  to  the  overall  system  by  constantly  changing  the  center 
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of  gravity  and  moments  of  the  combined  system . Because  the  moment  arms 
of  the  aerodynamic  forces  are  all  computed  about  this  center  and  because 
the  gravitational  force  is  applied  at  this  point , changing  these  quantities 
at  once  changes  all  the  torques  on  the  combined  system  from  what  they 
would  have  been  if  applied  to  the  projectile  body  alone.  In  addition,  since 
the  moments  of  inertia  are  constantly  used  in  the  rotational  equations  for 
the  total  system,  changing  them  will  at  once  also  change  all  these  equa- 
tions. Thus  when  the  combined  set  of  equations  is  solved  simultaneously 
by  the  computer  program,  the  internal  motion  is  described  by  Equation  18 
and  the  resulting  changes  in  the  center  of  gravity  and  moments  have  the 
proper  effect  on  the  total  projectile . 


All  of  this  can  be  summed  up  by  noting  that  the  relative  position  of 
the  total  center  of  gravity  in  the  missile  and  the  relative  velocity  of  this 
center  of  gravity  are  both  used  extensively  in  the  trajectory  program. 

Since  they  are  now  both  variable,  their  repeated  computation  must  be 
added  to  the  standard  analysis . The  following  sections  consider  in  more 
detail  the  modifications  relating  to  the  forces  and  moments  resulting  from 
the  movement  of  the  total  center  of  gravity  in  the  projectile,  in  order  to 
insure  that  this  mechanism  proceeds  correctly  in  the  trajectory  computation. 

Aerodynamic  Forces 


The  actual  aerodynamic  forces  acting  on  a projectile  depend 
(among  other  variables)  on  the  shape  of  the  projectile  body,  the  angle  of 
attack,  and  the  velocity  of  the  projectile.  The  expressions  for  the  various 
aerodynamic  forces  are  the  same  here  as  in  the  standard  trajectory  analy- 
sis and  they  will  not  be  repeated . But  the  velocity  they  depend  on  is  the 
velocity  of  the  projectile  body  or  skin  with  respect  to  the  air  V 
, ..  ' skin/air  ‘ 

n the  present  situation,  the  total  missile  is  not  rigid  due  to  the  movement 
of  the  internal  mass,  so  that  the  projectile  skin  and  the  total  center  of 
gravity  are  in  motion  with  respect  to  each  other.  Therefore,  the  velocity 

of  the  center  of  gravity  with  respect  to  the  air,  V which  was  former- 

cg/air, 

iy  used  in  the  aerodynamic  calculations,  is  no  longer  equal  to  V 
...  , skin/air 

which  should  more  properly  be  used,  so  that  it  becomes  necessary  to  com- 

PUtG  Vskin/airSeparately  for  the  combined  analysis.  However,  since  this 

difference  in  velocity  exists  only  in  the  component  along  th  H-system 
X-axis,  one  need  modify  this  component  only.  One  can  write  the  scalar 
equation,  assuming  XH-direction  components  only. 


V = v + V 

skin/ air  skin/cg  cg/air 


(19) 


15 


* 


Since  one  of  the  basic  functions  of  the  trajectory  program  is  to  compute 
the  motion  of  the  overall  center  of  gravity.  V . . can  be  calculated  using 

the  complete  system  of  equations  (including  Eq  18)  . The  first  term  on  the 
right,  is  calculated  using  Equation  12B,  noting  first  from  the 

figures  that  the  point  Q is  fixed  with  respect  to  the  projectile  skin  and  must 

move  with  it.  Since  the  distance  D is  measured  from  Q to  the  overall  cen- 

1 • 

ter  of  gravity,  measured  positive  along  the  + X -axis,  the  derivative  D 

H T 


as  given  by  Equation  12B  is  just  the  component  of  the  velocity  of  the  center 
of  gravity  with  respect  to  the  projectile  skin.  Thus  is  just  -D^, 

and  this  is  used  in  Equation  19  to  insure  that  the  proper  velocity  is  used 
in  computing  aerodynamic  forces. 


Generally,  these  forces  are  computed  by  multiplying  the  ex- 
pressions involving  velocity,  projectile  diameter,  and  air  density  by  pro- 
portional factors  called  aerodynamic  coefficients.  These  latter  quantities 
are  supplied  to  the  program  in  tables  dependent  upon  the  angle  of  attack 
and  Mach  number.  Since  the  Mach  number  used  for  this  purpose  should 
be  V , . . . (divided  by  the  speed  of  sound)  rather  that  V . . , the 
above  adjustment  in  velocity  given  by  Equation  19  is  again  necessary  to 
assure  that  proper  values  are  found  for  these  coefficients  themselves. 


Aerodynamic  Moments 

The  same  remarks  and  use  of  Equation  19  apply  to  the  compu- 
tation of  aerodynamic  moments,  which  similarly  depend  on  the  velocity  of 
the  skin  with  respect  to  the  air  and  involve  aerodynamic  coefficients  which 
do  aijo.  In  addition,  some  of  the  moment  terms  arise  from  the  fact  that  the 
aerodynamic  forces  act  through  a point  in  the  projectile  referred  to  as  the 
center  of  pressure.  Like  the  aerodynamic  coefficients,  these  centers  of 
pressure  are  also  found  from  tables,  themselves  functions  of  Mach  number. 
So  again,  the  remarks  above  relating  to  the  use  of  Equation  19  to  insure 
the  proper  values  are  relevant.  Finally,  since  these  moments  are  to  be 
calculated  about  the  overall  center  of  gravity,  which,  in  the  present  case, 
is  moving  in  relation  to  the  projectile  skin,  the  moment  arms  must  always 
be  computed  as  the  distance  between  the  center  of  pressure  and  the  current 
position  of  the  overall  center  of  gravity.  This  position  is  computed  using 
Equation  12A  first;  then  the  moment  arms  can  be  calculated  correctly. 


I 

1 


• i 
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Moments  of  Inertia 


i 


During  the  trajectory  calculation,  the  moments  of  inertia  of 
the  total  system  used  in  the  torque  equations  must  repeatedly  be  updated 
as  a consequence  of  the  motion  of  the  internal  mass.  Since  this  motion  is 
along  the  XH-axis  only,  the  moment  of  inertia  of  the  total  system  about 
this  axis,  denoted  by  I^iT)  , will  not  change,  but  the  moment  of  inertia 

about  the  YH-axis  Iy(T)  will  change.  Referring  to  the  fixed  moments  of 

inertia  of  the  internal  mass  and  body  alone  about  axes  parallel  to  the 
YH-axis  but  passing  through  their  own  centers  of  gravity  as  I^lm)  and 

1^  (B)  . respectively,  their  corresponding  moments  of  inertia  about  the 
YH-axis  (passing  through  the  total  center  of  gravity)  are  given  by 


I 


YH 


I 


YH 


(m)  = Iy  (m)  + m (D  - d)2 

(20A) 

(B)  = Iy(B)  + mB(DT  - Dfi)2 

(20B) 

Since  the  origin  of  the  H-system  is  always  located  at  the 
overall  center  of  gravity,  the  quantity  Iy(T)  , representing  the  moment 
of  inertia  of  the  total  projectile  about  a Y-axis  parallel  to  the  YH-axis  and 
through  its  own  center  of  gravity,  is  precisely  L_„(T)  . Also,  the  total 

projectile  is  merely  the  sum  of  the  internal  mass  and  the  projectile  body 
so  that  the  same  is  true  for  the  corresponding  moments  of  inertia.  To- 
gether these  imply  that 

VT)  = Iyh(T)  = IYH(m)  +Iyh(B)  (21) 

Substituting  Equations  20A  and  20B  into  Equation  21  yields  the  necessary 
result 


VT)  = IY(m)  + + m(DT  “ d)2  + mB(DT  ~ Db)2 

(22) 

which  is  used  to  compute  the  moment  of  inertia  of  the  total  system  using 
the  current  values  of  D^,  and  d.  It  should  be  remembered  that  because 

of  the  assumptions  of  symmetry,  l?  ( ) = I ( ) for  m,  B,  and  T,  so  that 
Equation  22  also  applies  with  Z substituted  for  Y. 


1 


' 
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Although  the  moments  are  updated  as  necessary  in  the  pro- 
gram using  Equation  22,  all  terms  involving  time  derivatives  of  moments 
of  inertia  have  been  neglected  in  this  study  (see  footnote,  page  9,  Ref  1)  , 
as  have  any  effects  of  rate  of  change  of  velocity  on  aerodynamic  forces. 

No  change  in  the  moment  I occurs,  so  that  1=0.  The  moment  Iv  (and 

also  I z)  varies  because  of  the  motion  of  the  internal  mass,  so  that  is 
actually  not  zero.  It  has  been  assumed  here  that  the  effect  of  the  deriva- 
tive term  I itself  is  small  compared  to  the  more  basic  effect  of  the  vari- 
ation of  I at  each  step  in  the  computation , which  is  accounted  for  by 

the  updating  using  Equation  22. 

The  present  analysis  has  been  limited  to  stages  of  trajectories 
without  fuel  burning,  so  that  no  mass  changes  occur.  Therefore  no  jet  re- 
action is  considered  (see  footnote,  page  7,  Ref  1)  and  no  changes  in  mo- 
ments due  to  mass  changes  occur. 

Application  to  the  Trajectory  Computer  Program 

From  the  standpoint  of  the  particular  trajectory  computer  program 
which  was  modified,  the  main  changes  consisted  of  adding  variables  repre- 
senting the  displacement  d and  its  derivative  d into  the  program  and  de- 
fining d by  Equation  18.  Connected  with  this  was  the  necessity  of  shifting 
the  total  center  of  gravity  (using  Eq  12A)  and  moments  of  inertia  (using 
Eq  22)  and  of  recomputing  the  velocity  of  the  skin  with  respect  to  the  air 
for  the  aerodynamics  (using  Eq  12)  in  the  innermost  loop  of  the  differen- 
tial equation  solving  code. 

Other  changes  involved  adding  new  inputs  describing  the  internal 
system,  converting  to  proper  program  internal  units,  initializing  variables 
and  minor  changes  to  the  output  to  provide  information  on  the  internal  sys- 
tem. The  alternate  set  of  equations  referring  to  the  vertical  (V)  coordinate 
system  (Ref  1)  had  previously  been  removed  so  that  consideration  of  this 
system  was  not  necessary. 

In  addition,  for  the  present  study  an  initial  "release  time"  was  im- 
plemented in  the  trajectory  program  so  that  one  could  specify  when  in  the 
course  of  the  trajectory  the  spring  would  be  released.  Until  this  time  is 
reached  in  the  flight,  the  internal  system  is  constrained  in  its  initial  posi- 
tion, as  specified  by  the  initial  condition  on  d.  When  this  time  is  reached, 
the  spring  is  in  effect  released  and  an  initial  condition  on  its  velocity  d 
can  be  applied . 


' 1 
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All  of  the  necessary  additions  and  modifications  based  on  this  anal- 
ysis have  been  implemented.  This  extended  an  existing  six-degree-of- 
freedom  missile  trajectory  program  at  Picatinny  Arsenal  to  accommodate  a 
missile  containing  an  internal  mass-spring  system  as  described  here. 

When  the  original  equations  plus  the  added  equation  of  motion  for  the  in- 
ternal mass  (Eq  18)  are  solved  simultaneously  (with  the  modifying  equa- 
tions described  in  the  discussion  of  the  combined  system)  , the  trajectory 
of  the  total  projectile  can  be  computed.  This  provides  both  the  trajectory 
of  the  combined  system  and  a solution  for  d and  d describing  the  position 
and  velocity  of  the  mass  within  the  projectile. 
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D„ 


D 

d 

d 


B 


eq 


m 


SYMBOLS 

Directed  distance  from  Q to  total  eg 
Directed  distance  from  Q to  body  or  skin  eg 
Directed  distance  from  Q to  eg  of  the  internal  mass 
Value  of  d at  spring  equilibrium  position 
Spring  constant 

Position  vector  of  internal  mass  in  H-system 
Position  vector  of  internal  mass  in  I-system 


m 


m 


B 


G 


B 


0} 

A 

V 

b 

gXH 


Vector  from  origin  of  I-system  to  origin  of  H-system 
Mass  of  internal  movable  object 
Mass  of  missile  body  or  Skin  alone 

Fraction  of  total  mass  attributable  to  missile  body  or  skin 
Force  vector 
Acceleration  vector 

Angular  velocity  vector  of  H-system  relative  to  I-system 
Velocity  vector 
Damping  coefficient 

Gravitational  acceleration  (iH  component  only) 

Unit  vectors  in  directions  of  H-system,  X,  Y and  Z-axes 


V 


skin/air 


Velocity  of  missile  skin  with  respect  to  air  (i  component  only) 
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V , • 

cg/air 

Vskin/cg 

Wm) 

iyh(b) 

iyh(t) 

Iy  (m) 

Iy(B) 

Iy(T) 


H 

• * # 

d,  d 

V 

XH 


dj  ( ) 

dt 


x 


Velocity  of  total  missile  eg  with  respect  to  air 
(ifJ  component  only) 

Velocity  of  missile  skin  with  respect  to  total  eg  (i^  component 
only) 

Moment  of  inertia  of  internal  mass  about  YH-axis 

Moment  of  inertia  of  missile  body  alone  about  YH-axis 

Moment  of  inertia  of  total  missile  about  YH-axis 

Moment  of  inertia  of  internal  mass  about  Y-axis  through  its 
own  center  of  gravity 

Moment  of  inertia  of  missile  body  alone  about  Y-axis  through 
its  own  center  of  gravity 

Moment  of  inertia  of  total  missile  about  Y-axis  through  its 
own  center  of  gravity 

NOTATION 


Refers  to  inertial  coordinate  system  which  may  be  considered 
to  have  its  origin  at  the  center  of  the  earth . 

Refers  to  horizontal  fixed-plane  coordinate  system , whose 
origin  is  at  the  center  of  gravity  of  the  total  missile. 


Dots  above  scaler  quantities  indicate  time  derivatives 

Component  of  vector  V in  direction  of  i (H-system  X-axis) 
This  component  may  also  be  indicated  by] 


XH 


Denotes  time  derivative  of  vector  realtive  to  I coordinate  system 


Denotes  vector  cross  product 
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